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Abstract. We do meson spectroscopy by studying the behavior of S'-matrix poles in the complex- 
energy plane, as a function of the coupling strength for ^Pq quark-pair creation. Thereto, a general 
formula for non-exotic hadron-hadron scattering involving arbitrary quark confinement is used, 
which can be applied to all flavors. We find two distinct types of poles, which we call confinement 
and continuum poles, respectively. Together, they suffice to understand the experimental meson 
spectrum. 



Introduction. From past studies on atomic spectra, we got so used to the terms "line 
spectrum" and "states" that we also apply them to mesonic spectra. Such a situation 
is depicted in Fig. la. The decay widths are thought of as higher-order effects, which 
indicate the life times of unstable states. This is indeed a fruitful strategy for atomic 
spectra. Starting from QED, one first determines the properties of the force which binds 
electrons to the nucleus, and next the corresponding bound-state spectra and states. 
Subsequently, one determines life times of states by considering the electromagnetic 
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FIGURE 1. (a): Spectra as we are used to. (b): The same, but restricted to one set of quantum numbers, 
(c): A modern spectrum of bound states below and resonances above threshold, in terms of complex poles 
of the scattering amplitude (Note that the units along the real and the imaginary axis are not the same in 
this figure). 

(EM) interaction of the photon field with the electrons. 

Seen from a different angle, the above strategy amounts to separating the interaction 
of electrons and the nucleus with the EM field into two distinct parts: the one binding 
electrons and the nucleus into an atom, and the one describing the scattering of photons 
off the atom. 

Here, we apply the same technique to strong interactions, but with the important 
difference that no manageable scheme exists to disentangle the two parts when setting 



out from QCD. Hence, they must be dealt with on the same footing. As a consequence, 
also the central resonance position and the resonance width must be described on the 
same footing, which can be achieved by referring to the complex pole position of the 
corresponding scattering amplitude. Such situation is depicted in Fig. Ic. 

Furthermore, the term "states" is not very well defined for resonances having large 
widths compared to their level separations. One may, of course, refer to the state at a 
certain energy, say the central resonance position. But it is not reasonable to assume that 
at energies say 100 MeV higher or lower, the state of the system is the same. Moreover, 
different resonances may be overlapping. 

Here, we shall be studying non-exotic meson-meson scattering. Through quark-pair 
creation, this implies multi-channel scattering, as well as intermediate states of confined 
systems of quarks and antiquarks. Hence, the advantage of complex spectra is apparent, 
since the pole structure of the full scattering matrix is the same for all channels, although 
the structure of peaks and zeros in the various channels can be very different. 

Meson- meson scattering amplitude. On rather general grounds [1, 2], it can be 
shown that a system of confined quarks and antiquarks coupled to free two-meson 
channels is well described by a partial- wave (£) inverse cotangent matrix of the form 
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Here, A parametrizes the intensity of quark-pair creation, p stands for the relativistic 
relative on-shell momentum in one of the two-meson channels, and E{p) is the total 
invariant mass of the system. Furthermore, En£ (n — 0, 1, . . .) represent the confinement 
spectrum, which is of course model-dependent, since it cannot be measured directly, 
while and .yV^i are convolution matrix integrals over quark-antiquark distributions, 
meson-meson scattering wave functions, and a transition potential. 

If we approximate the latter potential by a radial delta function, then we find for the 
ij-th matrix element of the partial-wave scattering amplitude 
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Here, jg represents the radial Bessel function, /i^^^ the radial Hankel function of the 
first kind, a the delta-shell radius, a/s the total invariant mass of the system, and jU, the 
reduced mass in the i-th two-meson channel. The coefficients r;(n) (vertex functions), 
which stem from the convolution integrals, are discussed later. 



Cross sections and poles. In Fig. 2 we show how cross sections following from 
formula (2) vary with increasing A, for S-wave isodoublet Kk scattering. In Fig. 2a 



the nonstrange-strange (ns) confinement spectrum is well visible for small X, whereas 
in Fig. 2c, for the model value of X, experiment is reproduced. 
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FIGURE 2. Cross section for 5- wave isodoublet Kk scattering. Left: For very small values of A, one 
observes the J'''^ = 0++ ns confinement spectrum. Middle: When X takes about half its model value, one 
notices some more structure for low invariant masses. Right: At the model's value of A, this structure 
is dominant and well in agreement with the experimental observations. The data are taken from Ref. [3] 

(open circles) and Ref. [4] (full circles). 

The pole structure for varying A is depicted in Fig. 3. In Fig. 3a we see how the 
pole stemming from the ground state of the confinement spectrum moves through the 
complex-energy (E = ^/s) plane when A changes. The pole movement shown in Fig. 3b 
is very different. For small A the pole disappears into the Kk continuum, which is the 
reason why it is not visible in Fig. 2a. 
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FIGURE 3. Pole positions for 5-wave isodoublet Kn scattering, as a function of the strength of quark- 
pair creation, parametrized by A (see Eqs . ( 1 ,2)) . (a) : Example of a confinement pole, which for A = lies 
on the real axis at total invariant mass Eq^, i.e., the ground state of the nsJ^ = 0+ confinement spectrum, 
(b): Example of a continuum pole, which for A = has an infinite negative imaginary part, i.e., the Kn 
continuum. 

The lowest-lying scattering-amplitude poles for DK in = O"*" and D*K in 7^ = 1 + 
are shown in Fig. 4. The two poles in Fig. 4a represent the £)*q(2317) resonance, which, 
in the absence of isospin-violating processes, comes out on the real axis, plus a broad 
resonance at about 2.8-2.9 GeV, undetected so far. 

The two poles shown in Fig. 4b represent the Z)^i(2460) and Dsi{2536) resonances, 
which both come out on the real axis. 



In Fig. 4a we depict two possible alternatives (one pair of solid lines and one pair of 
dashed lines). They result from very small variations in the delta-shell radius. Conse- 
quently, it is not at all trivial how to connect experimental poles with either confinement 
or continuum poles. 
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FIGURE 4. The lowest-lying poles for 5-wave = 0+ DK (a) and 7^ = r 
fat dots represent the pole positions for the model value of A . 

The two poles in Fig. 4b both stem from the ground state of the 7^ = 1 + charm-strange 
confinement spectrum, since there are two such states, i.e., and ^P\, which come out 
degenerate in our confinement spectrum. The linear combination 33% plus 67% 
decouples from D*K 5-wave scattering, hence becomes a bound state in the scattering 
continuum (see also Ref. [5]). On the other hand, the orthogonal linear combination of 
^P\ and ^P\ couples fully, and turns into a bound state below threshold. 

Threshold behavior. One may wonder how the lowest poles in Fig. 4 move along the 
real axis when X is increased. This is depicted in Fig. 5, where we continuously vary A 
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FIGURE 5. Movement of the D*()(2317) pole in 5-wave DK scattering near the DK threshold, as a 
function of A. Left: Movement in the p plane. Right; Movement in the E plane. We have marked the 
places where A equals 0.6, 0.7, 0.8, 0.9, 1.0, and 1.1. On the real axis in the E plane, the pole moves first 
towards threshold and then towards smaller real E for increasing A . 

from a small (0.60) towards a large value (1 .2). In the complex-momentum (p) plane, we 
observe that the pole turns purely imaginary somewhere in the lower half-plane. This is 
typical for S-wave scattering. For any higher angular momentum this happens at ;? = 0. 
The 5-matrix pole represents a bound-state when it lies on the positive imaginary axis. 



On the negative imaginary axis, it represents a virtual bound state. In this example, the 
pole lies on top of threshold for a value of X slightly larger than 0.7. 

The pole position in the E plane can be obtained from the pole position in the p plane 
through the relation 



Properties. The properties of the scattering amplitude following from the K matrix 
of Eq. (1) are the following. 

• Analytic in the total invariant mass. 

• Proper analytical and smooth behavior at all thresholds. 

• Unitary for real total invariant mass above and below thresholds. 

• Symmetric for any complex value of the total invariant mass. 

• Poles in the scattering matrix may be searched for in the same fashion as poles in 
any analytical function. 

• It contains an infinity of interfering resonances. 

• It also allows an arbitrary number (N) of coupled meson-meson channels. 

• For small A it exhibits the confinement spectrum 

• For the model's value of A, it reproduces the resonances and bound states of meson- 
meson scattering. 

• It contains kinematical Adler-type zeros [6] at almost the same invariant masses as 
the theoretical Adler zeros [7]. 

• It also reproduces automatically the low-lying scalar mesons. 

The parameters. The parameters A, a, and Ct) in formula (2) represent the probability 
of quark-pair creation, taken equal for up, down, strange and zero for the heavy flavors 
in the present version of the model, the average distance for quark-pair creation, and the 
average level spacing of the qq confinement spectrum, respectively. As expected, a is of 
the order of half a fermi, and 0) of the order of 200 MeV. In principle. A, a, and (a should 
follow from QCD. Furthermore, in our expression 



for the confinement spectrum, we must specify the constituent quark masses. We employ 
here the values m„ = 0.406 GeV, = 0.508 GeV, rric = 1 .562 GeV, and = 4.724 GeV 
fromRef. [1]. 

Moreover, flavor invariance demands that A^ and be scaled by the flavor mass, or 
the reduced quark mass in the case of different flavors [8, 9]. Consequently, for quark- 
pair creation, the average distance and probability decrease with increasing quark mass, 
as one naively would expect. 

Sum rules for vertex functions. The vertex functions, represented by r;(/i) in for- 
mula (2), where / indicates the channel index and n the radial excitation of the interme- 
diate resonance, contain the full structure of all quantum numbers, with no parameters 




(3) 



involved. The former can be determined from the recoupling matrix elements of the 
harmonic-oscillator wave functions for the effective quarks. The complete procedure is 
described in Refs. [10], assuming the ^Pq mechanism for quark-pair creation. 

Some of the resulting vertex functions are shown in Table 1, for a system with quan- 
tum numbers J^'^ — Each of the two scattering products (mesons) is characterized 

TABLE 1. Recoupling coefficients as a function of radial excitation n, for scalar 
(/ = 0, £ = 1,5= 1 , n) decay into two mesons. Notice that the recoupling coefficients 
for n = add up to one. This means that, in the harmonic-oscillator approach, 
there are no additional two-meson channels which can couple to the ground state 
of the confinement spectrum. For the higher radial excitations, the table is still very 
incomplete. 
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by its internal quantum numbers n£, J£, Lg, and Si (£ = 1,2), i.e., the internal radial ex- 
citation, total angular momentum (spin of the meson), orbital angular momentum of the 
qq system, and intrinsic qq spin, respectively. The relative motion of the two mesons 
is characterized by their relative orbital angular momentum L and total spin S, dictated 
by the /^^ quantum numbers. The relative radial excitation of the two-meson system is 
directly related to n (see Refs. [1 1] for details). 

When we determine the recoupling coefficients, we find that they decrease rapidly for 
higher radial excitations n (see Table 1), which implies that the higher terms in the sum 
over n in formula (2) are suppressed. 



5-wave scattering. In Fig. 2c we compare the result of formula (2) to the data of 
Refs. [3, 4]. We find a fair agreement for total invariant masses up to 1.6 GeV. However, 
we should bear in mind that the LASS data must have larger error bars for energies 
above 1.5 GeV than suggested in Ref. [4], since most data points fall well outside the 
Argand circle. Hence, for higher energies, the model should better not follow the data 
too precisely. 

Now, in order to have some idea about the performance of formula (2) for S'-wave 
7=1/2 Kk scattering, we argue that, as in our model there is only one non-trivial 
eigen-phase shift for the coupled Kn+Krj+Kri' system, we may compare the phase 
shifts of our model for Kr] and Kr]' to the experimental phase shifts for Kk. We do 
this comparison in Figs. 6 and 7, where, instead of the phase shifts, we plot the cross 
sections, assuming no inelasticity in either case. The latter assumption is, of course, 
a long shot. Nevertheless, we observe an extremely good agreement. In particular, for 
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FIGURE 6. 5-wave Kr] "cross section" (see text), as a function of total invariant mass, (a): From 
threshold up to 2.1 GeV. (b): Detail for lower energy, (c): Detail for higher energy. The data are taken 
from Ref. [3] (open circles) and Ref. [4] (full circles). 

Kt]' (Fig. 7) we become aware of a structure in the data at about 1.9 GeV, indicating the 
presence of a not- anticipated pole. This is something we would not have easily noticed 
from the data alone. 




FIGURE 7. 5- wave Krf' "cross section" (see text), as a function of total invariant mass. The data are 
taken from Ref. [3] (open circles) and Ref. [4] (full circles). 

When we inspect formula (2) for poles in the 5-wave isodoublet Kk scattering ampli- 
tude, then we find the pole structure as summarized in Table 2, i.e., five poles at energies 
up to about 2.2 GeV real part. The first pole, at 0.772 — 0.281/ GeV, describes the heavily 



disputed ^q(800) structure, whereas the second pole, at 1.52 — 0.097/ GeV, represents 
the well-established A^q(1430) resonance. 

TABLE 2. T-matrix poles for 5-wave Kn scattering, as obtained from Eq. (2). 



Pole (GeV) 


0.772-0.281; 


1.52-0.097? 


1.79-0.052? 


2.04-0.15? 


2.14-0.065? 


Origin 


continuum 


confinement 


confinement 


continuum 


confinement 



Our model is explicitly flavor-independent, meaning that the only flavor breaking in 
formula (2) stems from the effective quark masses which determine the ground state 
of the confinement spectrum (see Eq. 3), and from the masses of the mesons in the 
scattering channels. Consequently, titz scattering is not very different from Kn scattering 
in our model. We may expect then that each of the two flavor combinations that couple 
to isoscalar S-wave kk and KK scattering has a pole structure similar to the one in 
isodoublet Kn scattering, with the proviso that mixing in the 7 = case introduces 
an extra complication. The to-be-expected poles are given in Table 3, alongside the /o 

TABLE 3. The left-hand table shows the real parts of the poles in the ampUtude for nh (upper row) ss 

(lower row) 5-wave nn+KK coupled-channel scattering, as expected in analogy with Kn (see Table 2), 
whereas the right-hand table shows the observed /o states listed in the Tables of Particle Properties [12]. 
In the upper (lower) rows, we collect those resonances which we expect to be mainly nn (ss). 
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Structures reported in experiment. 

The often read comment that too many isoscalar states are observed [13], in order to 
justify the application of alternative quark, or even quarkless, configurations [14], is not 
confirmed here. 

Most probably, mesons are just mixtures [15] of quark-antiquark states [16], two- 
meson molecules [17], glueballs [18], tetraquarks [19], hexaquarks, hybrids [20], and 
so forth. Here, we have shown that the first two of the latter list of possible components 
are the most relevant ones. Moreover, a resonance is really a collection of states, all with 
different masses. Each of these states will have a different composition. 

In Table 4 we summarize the scalar poles which we obtained in the past [21, 22, 8]. 
The first four lines in the third column of Table 4 are refered to in the literature as 
the scalar-meson nonet [23], given by the isoscalars /o(600) and /o(980), the isotriplet 
ao(980), and the isodoublet pair 800) [24]. Here, they form part of the complex pole 
spectrum of the general scattering amplitude (1), where they appear as the lowest-lying 
continuum poles. The first four lines in the fourth column of Table 4 describe the nonet 
of scalar mesons that stem directly from the = 0"^ ground states of the confinement 
spectrum. 

Summary and discussion. In the foregoing, we have shown that the most relevant 

information on (meson) spectra amounts to the knowledge about the complex pole po- 
sitions of the scattering amplitudes. Furthermore, we presented a general form of am- 
plitudes for non-exotic hadron-hadron scattering with any possible flavor combination, 
based upon the most prominent properties of strong interactions, namely quark confine- 



TABLE 4. 5-wave scattering poles for various flavor/isospin 
combinations with 7^ = 0+. In the colunms "continuum", "ground 
state", and "excitation" we indicate the origin of the poles, i.e., ei- 
ther meson-meson continuum, or confinement spectrum (ground 
state or first radial excitation). 



channel 


q-q 


continuum 


ground state 


excitation 






GeV 


GeV 


GeV 




n-h 


0.47-0.21? 


1.36-0.13? 




Kn 


n-s 


0.77-0.28? 


1.52-0.10/ 


1.79-0.05? 


7]% 


n-h 


0.97 - 0.028/ 


1.45-0.13/ 




%% 


s-s 


0.99 - 0.020/ 


1.51-0.06/ 




Dn 


c-h 


2.14-0.16? 


2.58-0.12/ 




DK 


c-s 


2.33 


2.80-0.20? 




Bn 


n-h 


6.06-0.29/ 


5.46-0.03/ 


6.03-0.05/ 


BK 


s-h 


6.21-0.33? 


5.61 


6.05 - 0.03/ 


BD 


c-b 


7.12-0.43? 


6.64 


7.11-0.03? 



ment, quark-pair creation, and flavor invariance. The resulting pole spectrum may very 
well be compared to experiment. It is true that the predicted cross sections and phase 
shifts are not always in very accurate agreement with the data, as could hardly be ex- 
pected in view of the extremely wide scope of the model. In part this may also be due to 
the inaccuracy of some experimental numbers. Nevertheless, the model can be improved 
as well, as we suggest next. 

• At present, we determine the vertices in the framework of the ^Pq mechanism, 
assuming that all quark masses involved are the same. However, different quark 
masses can be dealt with, too. Computer programs already exist for several years, 
but have not yet been implemented. 

• The one-delta-shell approximation "unfortunately" works too well. Hence, al- 
though it is perfectly known how to implement more complicated transition po- 
tentials [25], and even some computer code is ready, the finishing touch could still 
take a while, unless some extra manpower becomes available. 

• The kinematics of meson pairs below threshold is most certainly not being dealt 
with in the most appropriate way by us. This is a subject which deserves to be 
studied in more detail. 
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